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Abstract 

We study a special anisotropic XYZ-model on a periodic chain of an odd length and 
conjecture exact expressions for certain components of the ground state eigenvectors. The 
results are written in terms of tau-functions associated with Picard's elliptic solutions of the 
Painleve VI equation. Connections with other problems related to the eight-vertex model 
are briefly discussed. 
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1 Introduction 



This is a sequel to our paper [1] devoted to connections of the eight-vertex model of statistical 
mechanics [2] with the theory of Painleve transcendents. Here we study a related case of the 

anisotropic XYZ- model on a periodic chain of an odd length, N = 2n + 1. 

(i) (i) (i) 

Let a x , Oy and a z , j = 1, . . . ,N, denote usual Pauli matrices acting at the j-th site of 
the chain. Consider a particular XYZ-Hamiltonian 

N 
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1 N 

hxyz = ~E( J * + jy + (i) 

where the coefficients 

_ 2(1 + C) , _ 2(1-C) r _ C^l 
Jx ~ C 2 + 3 ' v ~ + 3 ' z ~ C 2 + 3' 1 j 

are specific rational functions of a single parameter which satisfy the relation 

Jx Jy + Jy Jz + Jz Jx = . (3) 

Baxter [3] proved that for an infinitely large chain, N — ► oo, the ground state eigenvalue of ([T|) 
in this case has a very simple form 

£q(Q Jx + Jy + Jz _ 1 , , 

iV 2 2' U 

Later on it was conjectured [4] that this expression is exact for all finite odd values of N. 

The Hamiltonian ([1]) commutes with the transfer matrix of the eight-vertex model (8V- 
model), where the Boltzmann weights a, b, c, d (we use standard notations of [2], see Sect. [2] 
below for further details) are constrained as 

(a 2 + ab) {b 2 + ab) = (c 2 + ab) {d 2 + ab) , (5) 

and the variable £ in ([2D is given by 

cd (a — b + c — d)(a — b — c + d) 

^ = d> , 7= (a + b + c + d)(a + b-c-d)' ^ ' 

Additional variable 7, introduced for later convenience, is connected to ( by a simple self- 
reciprocal rational substitution 

,7 + 3 C + 3 

c = — ' 7= c^T- (7) 

The spectrum of (HD possesses an S3 symmetry group with respect to permutations of the 
constants J x , J y and J z . Indeed, any such permutation can be compensated by a linear transfor- 
mation which acts on the eigenvectors and does not affect the spectrum. For the parametrization 
([2]) this group is generated by two substitutions of the variable 



S X y . £ ► £ ',' J X <— * Jyi Jz > J%i 

s ■ <T - => J ~ J J -J (8) 

°xz • S £_i ' J i I'zi "2/ ' J j/ • 
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The largest eigenvalue of the transfer matrix, corresponding to ([!]), also has a remarkably 
simple conjectured form [4] 

A = (a + b) N , N = 2n + 1, (9) 

which is expected to hold for finite chains. 

In [5] we studied Baxter's famous TQ-equation for this simple eigenvalue and found 
corresponding eigenvalues of the Q-operator. With an appropriate normalization they can be 
expressed through certain polynomials 



V n (x,z) = Y j rt\z)x\ z = T\ n = 0,l,2,... , (10) 



k=0 

of the variable x, defined by the following quadratic equation (see also ([23]) ). 

7 \2 16 (a - b) 2 cd 



y/xJ (c + d) 2 (a + b + c + d)(a + b-c-d)' 



11) 



The coefficients (z), i = 0, . . . ,n, appearing in (jlOH . are polynomials in the variable z = r y~ 2 
with positive integer coefficients. Detailed definitions of the polynomials (|10p are presented 
in Sect. [2j Here we want to illustrate their connection to the Painleve VI equation. This 
connection manifests itself in some specific properties of the coefficients rj (2). In particular, 
let s n {z) = rn(z) be a coefficient in front of the leading power of x in (|10p . In [5] we conjectured 
the following recurrence relation 

2z{z - 1)(9j8 - l) 2 d 2 z logs n (z) + 2(3z - l) 2 (9z - l)d z logs n {z)+ 

+8 ( 2n + 1 )2 s n+l(^n-l(^) _ [4(3n + 1)(3n + 2) + (£)z _ 1)n(5n + 3)] = Qj ( 1 2 ' 

where sq(z) = s\(z) = 1, which uniquely determines the polynomials s n (z), for all n£Z. Later 
on we proved [1] that Eq. (|12p exactly coincides with the recurrence relation for the tau- functions 
associated with special elliptic solutions of the Painleve VI equation. In this letter we extend 
these connections to study ground state eigenvectors of Hxyz> corresponding to the eigenvalue 

©• 

For odd N all eigenvalues of (pQ) are double degenerate. Thus, there are two ground state 
eigenvectors 

Hxyz*± = ^o^±, 5^± = ±^±, Ky± = y T , (13) 

where 

S = o-M ® a f) ® . . . ® CT f ) 5 n = a (i) 8 ^(2) g, . . . 3 CT (7V) ^ (M) 

and 

[Hxyz,cS] = [Hxyz, ft] = 0, nS = (-l) N SK. (15) 

Due to the spin reversal symmetry, generated by the operator 1Z, it is enough to consider one 
of these vectors. For definiteness, consider the vector , I / _. Omitting the suffix "— ", we denote 
its components as ^i lt i 2 ,...,i N , where ii,i2, ■ ■ ■ , *jv £ {0, 1} and assume an orthonormal basis \i), 
i = 0, 1, for each spin 

ff*|0} = +|0), a z |l)=-|l). (16) 

The ground state eigenvectors are translationally invariant and possess a left-right reflection 
symmetry. Taking this into account we will give only one non-vanishing representative from 
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each symmetry class. Note that for non- vanishing components ^i u i 2 ,...,i N of \P_ the number of 
"down-spins" in the set {ii,i2, ■ ■ ■ ,in} is odd, while for vanishing components it is even, 

*M 2 ,..,iiv = 0, if h+i 2 + --- + iN = (mod 2). (17) 

The fact that both the coefficients ([2]) and the eigenvalue (jl]) are rational functions in £ with 
integer coefficients implies that with a suitable normalization one can make all components 
of the eigenvector , 4?_ to be polynomials in £ with integer coefficients [6] (such that there are 
no polynomial factors common for all components). This choice is unique up to a numerical 
normalization. The later is fixed by the requirements 



,0 



= 1, for odd n; 01 



= 1, for even n. (18) 



Note that in the case £ = the Hamiltonian (TTJ) reduces to that of the XXZ-model with the 
parameter A = —1/2. From this point of view the normalization (|18j) is identical to that used 
in [7] j where this particular XXZ-model was studied. 

We have calculated all components of the eigenvectors directly from the definition (|13p for 
N < 17 (and some particular components for N < 25) and made several interesting observations 
which we formulate as conjectures valid for all N = In + 1. As an example we present here 

Conjecture 1. The norm of the eigenvector is given by 

l*-| 2 = E K,i,,.,i, = (VS)" C (n+1) Sn(C 2 ) S-n-l(C 2 ) , (19) 

ii,i2,---,ijve{0,l} 

where s n (C~ 2 ), n £ Z, are defined by the recurrence relation with z = £~ 2 and sq(z) = 
si(z) = 1. 

Other conjectures on the properties of eigenvectors require additional notations; they are 
presented in Sect. [3l Basic definitions for the 8V-model, a brief review of some of our previous 
results [5] and some new results on the eigenvalues of Baxter's Q-operators are given in Sect. [2j 
In Conclusion we discuss some unresolved questions and connections of our results to other 
problems related with the eight- vertex model, in particular, to the eight- vertex solid-on-solid 
model [8] with the domain wall boundary condition [9] and the three-coloring problem [10]. 



2 The eight-vertex model and TQ-equation 

2.1 Basic definitions and notations 

We consider the eight-vertex model on the iV-column square lattice with the periodic (cylindri- 
cal) boundary conditions and assume that N is an odd integer N = 2n + 1. Following [2] we 
parametrize the Boltzmann weights a, b, c, d of the model aaj 

a = p #4(277 1 q 2 ) -di(u — 77 1 q 2 ) "9\(u + 77 1 q 2 ), 



1 We use the notation of [11] for theta-functions $fc(«|q), k — 1,2,3,4, of the periods n and ttt, q = e 
Imr > 0. The theta-functions H(v), Q(v) of the nome qs used in [2] are given by 

qs = q 2 , h(«) = <m^-| q 2 ), e(«) = 04(^-1 q 2 ), 

where Kb(&) is the complete elliptic integral of the first kind with the elliptic modulus k = - !?2(0|qB)/i?i(0|q_B) 
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b = p #4(277 1 q 2 ) #i(u - 77 1 q 2 ) # 4 (u + 77 1 q 2 ), 

c = p #1(277 1 q 2 ) # 4 (u - 77 1 q 2 ) # 4 (u + 77 1 q 2 ), (20) 
d = p #1(277 1 q 2 ) #i(u - 77 1 q 2 ) #i(u + 77 1 q 2 ), 

and fix the normalization factor, 

p = 2 #2(0 I q)- 1 # 4 (0|q 2 )- 1 . (21) 

With this parametrization the constraint © is equivalent to the condition 

77 = vr/3 , (22) 

which will be always assumed throughout this paper. This still leaves two arbitrary parameters: 
the (spectral) parameter u and the elliptic nome q = e l7TT , Imr > 0. The variables £, 7 and x, 
defined in ©, and (fTTj) . can be written as 



r#i(f |q 2 )" 


to 


\M% lq 1/2 )" 


to 


r# 3 (f Iq 1 / 2 )! 




> 7 = - 


, x = 7 




k(f |q 2 )_ 


U 2 (f Iq 1 / 2 ). 


L# 4 (||q 1 / 2 )J 



z = 7 ' 2 • (23) 



Note that the last expression for x determines our choice of a particular root of the quadratic 
equation (fTTI) . 

2.2 The TQ-equation 

Any eigenvalue, T(u), of the row-to-row transfer matrhj^l of the 8V-model satisfies Baxter's 
famous TQ-equation [2], 

T(u) Q(u) = <j)(u - 77) Q(u + 277) + 4>(u + 77) Q(u - 2 V ), (24) 

where, with an account of (|2ip . 

4>{u)= #f(« I q). (25) 
With the parametrization (|20p . (|2ip the eigenvalue ([9]) takes the form 

T(u) = {a + b) N = <p(u), 77 = vr/3, iV = 2n + 1. (26) 

Equation (f24"j) for this eigenvalue, has been studied in [5]. It has two different solutions [12,13], 
denoted Q±{u) = Q±(u, q, n), which are entire functions of the variable u and obey the following 
periodicity conditions [2,14], 

Q±(u + tt) = ±(-l) n Q±(u), Q±(u + 7tt) = q"^ 2 e ~ iNu Q T (u), Q±{-u) = Q±(u). (27) 

The above requirements uniquely determine Q±(u) to within Tt-independent normalization fac- 
tors. The solutions Q±{u) satisfy the quantum Wronskian relation [12,15] 

Q+(u + rj) Q_(u - rj) - Q+(u - r]) Q.(u + rj) = 2i<p{u) W(q, n), (28) 

where W(q,n) is a function of q and n only (The fact that W(q,n) does not depend on the 
variable u follows from (|24p and (|27p ). Note that, taking into account the periodicity (|27p . one 
can bring Eq. (|24|) to the form 

<p(u) Q{u) + cf)(u + 2vr/3) Q(u + 2vr/3) + <f>(u + 4vr/3) Q(u + 4vr/3) = . (29) 
2 We use exactly the same definition of the transfer matrix as in [2]. 
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Below it will be more convenient to use the combinations 

Qi(«) = (Q+(u) + Q-(u))/2, Q 2 {u) = (Q+(u) - Q_(«))/2 , (30) 
which are simply related by the periodicity relation 

Q$(u + 7r) = (-l) n Q$(tO. (31) 

Bearing this in mind we will only quote results for Qi(u), writing it as Q^(v) to indicate the 
n-dependence. Introduce new functions V n {u) instead of Q^fu), 

Qf\u) = jV(q, n) | q 1 / 2 ) tf 4 2 >/2 | q 1 / 2 ) V n (u), (32) 

where jV(q,n) is an arbitrary normalization factor. The analytic properties of V n (u) are de- 
termined by the periodicity relations ([2"T|) and the fact that the eigenvalues Q^\{u) are entire 
functions of the variable u. A simple analysis shows that V n {u) is an even doubly periodic 
function of u, 

Vniu) =V n (u + 2vr) = V n {u + ttt), V n {u) = V n {-u) , (33) 

with all its poles^l (of the order 2n and lower) located at the point u = ttt/2. Every such 
function is an n-th degree polynomial in the variable x, given by (|23|) (see § 20.51 of ref. [11]). 
The coefficients in these polynomials will, of course, depend on the elliptic nome q. Let us now 
change independent variables from u and q to the variables x and z = r y~ 2 , defined in ([23"j) . and 
(with a slight abuse of notations) write V n (u) as 



V n (x,z)=Y,r < k l \z)x k . (34) 



k=0 

The TQ-equation f|29|) can be re- written in terms of the polynomials V n (x,z). For a fixed 
value of the nome q, the variable x in (|23|) is a function of u, so we can write it as x = x(u). 
Introduce two new variables 

x± = x(u± f) = 7 2 /x(n=F )• (35) 

They satisfy the relations 

(x - l) 2 _ 2z(x 2 z + l)-x(z 2 + Az-l) 

%+X- — , , w i a; + + x - — , -, \o j (3DJ 

(x z — l) z z (x z — lj z 

which can be easily solved for x± in terms of x and z. The resulting expressions involve a square 
root from a third order polynomial in x. It is convenient to define 

1 x(z-l)[(2x-3)z + l] x±-l 
J± 2 2z(x_-x+)(xz-l) 2 ' P± (l-zx ± )x' 1 ' 

With all these new notations the TQ-equation ()29j) can now be transformed to its algebraic form, 

V n (x, z) = p + f 2n+1 Vniz^xZ 1 , z) + p_ /f +1 ^(z-^; 1 , z) . (38) 

Substituting (|34p into the last equation, and expanding it near the point x = 0, one imme- 
diately obtains a simple relation, 

r ( ®\z) = V n (0, z) = 4~ n z- 1 (z + n(3z - 1)) V n { Z -\z) - A~ n z~ 2 (z - 1) dVn ^ z) , (39) 

OX x=z~ 1 

quoted here for future references. 



3 An apparent pole at u — ty + ttt/2 cancels out because Q^iu) vanishes at this point as a consequence of the 
periodicity conditions (J27J. 
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2.3 Polynomials V n (x, z) 

Let us now substitute polynomials (|34p into the TQ-equation (|38p . Excluding x+ and :c_ with 
the help of (|36p . one can readily see that the RHS of (|38|) is a rational function of s (indeed, it 
is a symmetric function of x+ and x_ and, therefore, can be expressed through two elementary 
combinations (I36p ), Writing Eq. (|38p as a polynomial in x and equating its coefficients to zero 
one obtains an (overdetermined) system of homogeneous linear equations for n + 1 unknowns 
r^\z), r^\z), . . . ,r!^\z). All elements of the coefficient matrix for this system are rational 
functions of the variable z with integer coefficients. This means that with a suitable normal- 
ization all r^\z), k = 0, 1, . . . ,n, can be made polynomials in z with integer coefficients (such 
that there are no polynomial factors common for all r^\z)). Thus, V n (x , z) are two- variable 
polynomials in x and z with integer coefficients. The first few of them are given in (|46p and 
Appendix A below. 

Originally, we have calculated [5] these polynomials for n < 10 by directly solving Eqs. (|29p 
and (|38p by a combination of analytical and numerical techniques. For larger n this did not 
appear be to practical due to complexity of intermediate expressions. Subsequently, in the same 
paper [5], we found a more efficient method for the calculation of V n {x, z), based on the partial 
differential equation (|4*4"j) . discussed below. We have observed that all coefficients of V n (x,z) 
are, in fact, positive integers for all n < 100 and suggested that these coefficients might have a 
combinatorial interpretation (which is yet to be found). 

Below we summarize all important properties of T , n (x,z) discovered in our previous works 
[1,5]. 

Conjecture A ( [1,5]). 

(a) The degrees of the polynomials r^ n \z), i = 0, . . . , n, appearing as coefficients in the expansion 
(34\ ), are given by 

deg[rj; n) (*)] < [n(n - l)/4 + k/2\ , (40) 
where [x\ denotes the largest integer not exceeding x. 

(b) if the normalization ofP n (x,z) is fixed by the requirement 

ri»>(0) = 1 , (41) 

then all polynomials (z), k = 0,1,..., n, have positive integers coefficients in their 
expansions in powers of z. 

The normalization (|4ip will be implicitly assumed throughout the rest of the paper. The 
most important property of the polynomials V n {x^z) is that they satisfy a remarkable linear 
partial differential equation. This equation can be written in different forms, depending on the 
choice of independent variables and unknown function. First consider the case of the original 
variables u and q of the 8V-model. Introduce the functions 

o2n+l/ I \ 

*±(«,q,n) = 1 \ ' 3 q ; Q±(n,q,n), (42) 

where Q±(u, q, n) are eigenvalues of the Q-operators, defined in Sect l2.21 The analytic properties 
of 3>±(it, q,n) in the variable u are determined by (|27p . 
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Conjecture B ( [5]). The functions $±(u, q,n), defined by satisfy the non- stationary 

Schrddinger equation 

d d 2 
6q—$(u,q,n) = {- — + 9n(n + l) P(3u \ q 3 ) + c(q, q, n). (43) 

Here the modular parameter r plays the role of the (imaginary) time and the time-dependent 
potential is defined through the elliptic Weierstrass p-function [11] (our function p(v \ e t7re ) has 
the periods it and ire). The constant c(q,n) appearing in (|43p is totally controlled by the 
normalization of Q±(u) and can be explicitly determined once this normalization is fixed (see 
Eqs.(37) and (38) in [5]). Equation (j4"3|) is obviously related to the Lame differential equation 
and could be naturally called the "non-stationary Lame equation". This equation arises in 
various contexts [16, 17] which are not immediately related to this paper. 

The differential equation (|43h can be equivalently rewritten in an algebraic form for the 
polynomials V n {x, z). 

{A(x,z)d 2 x +B n (x,z)d x + C n (x,z) +T(x,z)d z }v n (x,z)=0, (44) 

where 

A(x, z) = 2x(l + x — 3xz + x 2 z)(x + 4z — 6xz — 3xz 2 + 4x 2 z 2 ), 
B n (x,z) = 4(1 + x - 3xz + x 2 z)(x + 3z - 7xz + 3x 2 z 2 )+ 

+2nx(l - Uz + 2lz 2 - 8x 3 z 3 + 3x 2 z(3z 2 + 6z - 1) - x(l -9z + 23z 2 + 9z 3 )), 

(45) 

C n (x,z) = n[z(9z - 5) +x 2 z(3z 2 + llz - 2) + x(9z 3 - 38z 2 + 19 z - 2)- 

-4x 3 z 3 + nz(l -9z- x(9z 2 - 36z + 3) + x 2 (3z 2 - 3\z + 4) + 8x 3 z 2 )}, 
T(x,z) = -2z{\- z)(l-9z){l + x-3xz + x 2 z). 

It is fairly easy to prove [5] that the differential equation (|4"3|) , restricted to a class of functions 
<J?(«, q,n) with suitable analytic properties in the variable u, implies the functional equation 
(|29p . The non-trivial part of the Conjecture B is the fact of existence of solutions of (|43p with 
these analytic properties. For Eq. (|44p this translates into a question of existence of solutions, 
which are polynomials in the variable x. 

Equation (|44p is extremely useful for finding polynomial solutions, even though the coeffi- 
cients therein look very complicated. The first polynomials 'P n (x,z) read 

V (x,z) = l, Vi(x,z)=x + 3, V 2 (x,z) = x 2 (l + z) + 5x(l + 3z) + 10, 

, (46) 

V 3 (x, z) = x 3 (l + 3z + Az 2 ) + 7x 2 (l + 5z + 18z 2 ) + 7x(3 + 19z + 18z 2 ) + 35 + 21z, 

the next one is given in Appendix A. The constant term and leading coefficient in these poly- 
nomials (with respect to the variable x) are determined by the following 

Conjecture C ( [1,5]). The coefficients for the lowest and highest powers of x in V n (x,z), 
corresponding to k = and k = n in ft34\ ), read 

= ^\z) = r n (z, -1/3), s n (z) = rW(z) = r n+1 (z, 1/6) , (47) 
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where the functions T n (z, £) (for each fixed value of the their second argument are determined 
by the recurrence relation 

2z(z - l)(9z - l) 2 [logr„(z)]'; + 2(3z - l) 2 (9z - l)[\og T n (z)]' z + 

l T n+ l(z)T n -l(z) 



2n - 4£ - - 
S 3 



T n( z ) 

-[12(3n - 6£ - l)(n - 2£) + (9z - l)(n - l)(5n - 12f)] = 0, (48) 
with the initial condition 

r (z,O = l, n(z,0 = -4^ + 5/3 . (49) 
T/ie functions T n (z, £) are polynomials in z for all n = 0, 1, 2, . . . , oo. 

As explained in [5], the partial differential equation leads to a descending recurrence 
relations for the coefficients in (|34p . in the sense that each coefficient r^, n ^(z) with k < n can 
be recursively calculated in terms of (2), with m = fc + 1, ...,n and, therefore, can be 
eventually expressed through the coefficient r^\z) of the leading power of x. Conditions that 
this procedure truncates (and thus defines a polynomial, but not an infinite series in negative 
powers of x) completely determine the starting leading coefficient as a function of z. The above 
conjecture implies that these truncation conditions are equivalent to the recurrence relation (]12p 
which is a particular case of (|49|) for £ = 1/6. Similar reasonings apply to the coefficient Tq(z) 
in (|34p (the constant term with respect to the variable x). 

Note, that Eq. (|48p exactly coincides [1] with the recurrence relation for the tau- functions 
associated with special elliptic solutions of the Painleve VI equation. 

2.4 Quantum Wronskian 

We conclude this section with a short analysis of the algebraic form 

V n {x + ,z)V n (z~ 1 xZ 1 ,z) + ; — ±- -V n {x-, z)V n {z~ 1 %+ 1 , z) 

( ,-i) V J Wn(z) 

(50) 

of the quantum Wronskian relation (I28p . Here \N n (z) is related to W(q,n) in (|28p . 

W(q,n) = (-l) n z[2i9i(7r/3|q)] 2ri+1 AA 2 (q,n)W n (2;) . (51) 

Equation (|50p is an algebraic identity valid for arbitrary values of x. Expanding this identity 
around x = z -1 , one obtains 

\N n (z) = -s n (z)V n (z- 1 ,z), n>0, (52) 

where s n (z) is defined by (fT2"|) (it coincides with the Painleve VI tau- function s n (z) = T n+ i(z, 1/6) 
defined by (|4"8"j) for ^ = 1/6) . Similarly, expanding ([50]) around x = and using ([52]) , one obtains 

^„(1,«) =4 n s n (z), n>0. (53) 

Interestingly the quantity V n (z~ 1 , z), entering (]52p is also determined by the Painleve VI recur- 
rence relation (1481). 



1 — X + X-(l — xz) ' ' 1 — x + x+(l — xz) 

1 ( z(xz — l) 2 {x + — x_) 2 
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Conjecture D. The value V n {z 1 ,z) is determined by recurrence relation J^#| j with £ = 2/3, 

V n {z-\z) = (-4/3)" z~ n T n+2 (z, 2/3), n > 0. (54) 

Combining the above formulae one obtains the following expression for the quantum Wron- 
skian, 

\N n {z) = -(-4/3)" z- n r n+1 {z, 1/6) r n+2 (z, 2/3), n > . (55) 

3 Eigenvector results 

To formulate our results for the eigenvectors (|13p we need to define an additional set of poly- 
nomials p n (y) and q n (y), n £ Z. In principle, these polynomials can be defined by yet another 
recurrence relation of the Painleve VI type (though more complicated than (|48p ) which will be 
presented elsewhere. For our purposes here it is much simpler to define these new polynomials 
p n (y) and q n {y) as subfactors of already introduced polynomials s n (z). We will do this by means 
of the Conjecture E, given below. Remind that s n (z), n £ Z are defined by Eq. (fl~2j) with the 
initial conditions so(z) = si(z) = 1. 

Conjecture E. 

(a) The polynomials S2fc+i(y 2 ) factorize over the integers, 

S2k+i(y 2 ) = s 2 k+i{o)Pk(y)Pk(-y), Pk(fy = l, k e z, (56) 

where Pk(y) are polynomials in y with integer coefficients, degp^iy) = k(k + 1), such that 
p' k (0) > 0, k > 1 and p' k (0) < 0, k < —2, where p' k (y) = dpk{y)/dy denotes the derivative in 
y. Note that p-i(y) = p (y) = 1. 

(b) the polynomials Pk(y) possess the symmetry 

/ l + 3y \*(fe+i) ( l_-y \ 
Pk(y) = { 2 ) Pk \ l + 3y )' ' ^ ' 

(c) polynomials S2fc(y 2 ) factorize over the integers, 

S2k{y 2 ) = c k (1 + 3y) fc ( fc+1 ) (f^-) Qk-i(y), k e Z, (58) 

where qt{y) are polynomials in y with integer coefficients, degqk(y) = k(k + 1), Q'fc(O) = 1 
and 

c k = 2-^ k+2 \ k>0; c k = 2- fc2 (2/3) 2fc+1 , k<0. (59) 

(d) polynomials q k {y) possess the symmetry 

/l + 3y\fc(fe+i) ( y~l \ , c „ 

g*(y) = ( 2 j gfc (, 1 + 3j/ J' (60) 

A few first polynomials s n (z), p n {y) and q n (y) is listed in Appendix A. Let us mention one 
simple, but important corollary of the above conjecture. The LHS of (|58p is an even functions 
of the variable y. Combining this fact with the symmetry relation (|57p . one immediately deduce 
that q k {y) is also an even function, 

Qk{-y) = qk(y), k e z . (61) 
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In Conjecture 1, given in the Introduction, we have stated an explicit expression (|19|) for the 
norm |\£_| 2 of the eigenvector f _ as a function of the parameter £ entering the Hamiltonian (pQ). 
Note, that using the factorization and symmetry properties (|56p -(i6"T1). one can show that the 
rescaled norm 

N n (() = (C 2 + 3)" n(n+1)/2 |^_| 2 (62) 

is invariant with respect to the full S3 symmetry group generated by the substitutions ([8]), 
which is a well expected result. Further, is easy to see, that modulo a trivial numerical factor, 
the expression for the norm remains unchanged upon the replacement n — * — n — 1, which 
corresponds to a negation of the length of the chain, N — ► —N. In other words, the norm is an 
even function of the length of the chain. It would be interesting to understand a reason of this 
phenomenon. 

We are now ready to present further conjectures on the properties of the eigenvectors. 
Conjecture 2. The component of the eigenvector \E'_ with one spin down is given by 

V>o...ooi = ^C {n - 1)/2 UC 2 ), N = 2n + 1, (63) 

where s n (z) is defined by ffl 7 ^ . 

A few first polynomials s n {z) and s n (z) is listed in Appendix A. 
Conjecture 3. The component of the eigenvector with all spins down is given by 

V>U...U = C [n+1)/2 s n (C 2 ), N = 2n + 1 , (64) 

where s n (z) is defined by 

It is interesting to note that to within a simple power of £ the above two components of the 
eigenvector precisely coincide with the constant term and leading coefficients of the polynomial 
V n (x,z), which is simply connected (|32p with the corresponding eigenvalue of the Q-operator. 
We believe that this fact certainly deserves further studies. 

Finally, consider components of with alternating (up and down) spins in the chain, 

A n (C) = *ooioi...oi, for odd n; A n (() = *oioi...on for even n . (65) 

In the case £ = these are largest components of the eigenvector. 
Conjecture 4. The components of with alternative spins are given by 

A 2k (C) = 2 M2-*) (3 + c) fc ( fc -V (fe - 1 Wife) gjk-xCr 1 ) 

(66) 

A 2k+ i(C) = 2- fc2 (3 + C) fc ( fc+1 )C fc(fc - 1) P fc (^) Qk-iiC 1 ) 

A few polynomials p n (y), Qn(y) an d A n (Q is listed in Appendix A. As noted before, the 
case C = corresponds to the XXZ- model with A = —1/2. It is known [7, 18-20], that in this 
particular case, the values of the components (|65p . normalized by (|18p . coincide with the number 
of alternating sign matrices 

k=0 y ' 
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calculated in [21]. Using this result in (|66|) one easily obtains for n > 0, 



Vn 



K (2\ n{ - n+1) ^ (2fc)!(6fc + l) 



1 ' X JL (4fc)!(4fc + 1)! 



n(n+1) , x,, _ n _! A (2fc + l)!(6fc + 4)! 
O V(C )k=o - 2 11(^ + 2)1(4* + 3)! ■ (69) 

Apparently one can derive these expressions directly from the definitions of the polynomials 
Pn{y) and q n (y), given in Conjecture E, however, we postpone this to a future publication. 

Finally, mention one amusing observation connected with the expressions (|66p . It is not 
difficult to analytically derive an asymptotic expansion A n asymp \Q which correctly reproduce 
first terms of the expansion of A n ((^) for small £ up to the order 0(£ 2n ), 

A n (() = At ymp) (0 + G(C 2n ), C - . (70) 

Analytically continuing this asymptotic expansion to n = 0, 

A (asymp) (^ = 1 _ ^2 _ 3 ^4 _ 15 ^6 _ gg ^8 _ 534 ^10 _ 3473 ^12 _ _ _ _ ^ 

n=0 

and plugging its coefficients into Sloane's integer sequences database (in a search for a discovery) 
we found that they only "slightly" mismatched numbers of lattice animals made of n three- 
dimensional cubes [22], which are 1,3,15,86,534,3481,.... Of course, it would be extremely 
weird if they matched. 



4 Conclusion and outlook 

In this paper we have demonstrated that a particular anisotropic XYZ-model, defined by ((T|) 
and d2|), is deeply related with the theory of Painleve VI equation. We have proposed exact 
expressions for the norm (Conjecture 1) and certain components of the ground state eigenvectors 
(Conjectures 2,3,4). The results are expressed in terms of the tau- functions associated with the 
special elliptic solutions of the Painleve VI equation [1]. 

In this connection, it is useful to mention other celebrated appearances of Painleve transcen- 
dents in mathematical physics. The most prominent examples include the two-dimensional Ising 
model [23], the problem of isomonodromic deformations of the second order differential equa- 
tions [24] and the field theory approach to dilute self-avoiding polymers on a cylinder [25-29]. 
The latter problem is connected with the massive sine-Gordon model at the supersymmetric 
point (where the ground state energy vanishes exactly due to supersymmetry). Our previous 
work [5] grew up from attempts to develop an alternative approach to this polymer problem 
based on the lattice theory. It turns out that all non-trivial information about dilute polymer 
loops is contained in the ground state eigenvalues [5] of the Q-operator for the 8V-model on a 
periodic chain of an odd length, connected with the special XYZ-model, considered in this paper. 
In [5] we have found that these eigenvalues can be uniquely determined as certain polynomial 
solutions T , n (x, z) of the partial differential equation (I44D (remind that the variable x and z are 
connected to the original spectral parameter u and the elliptic nome q, respectively, see (|23p ). 
So far we have not ultimately understood the role of this equation in the Painleve VI theory, but 
there is no doubt that there are profound connections. For example, one-variable specialization 
of V n {x, z) at particular values of x (which remain polynomials in the variable z) are connected 
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with the tau-functions associated with the Picard solutions of the Painleve VI equation (see 
Eqs. ([53|) and (JM|))- The same property is enjoyed also by the coefficients in the expansion of 
V n (x, z) in powers of x. Note, that namely these coefficients provide "construction materials" in 
the expression for the ground state eigenvectors (see Eqs, (|19p . (|56 p - (|66p ). Most of our results are 
conjectures and it is, of course, desirable to obtain their proofs. Another outstanding problem 
is an algebraic construction of the Q-matrix. As noted in [30], the method used in [2] for the 
construction of the Q-matrix cannot be executed in its full strength for ij = it/ 3, since some 
axillary Q-matrix, Q/j(u), in [2] is not invertible in the full 2 Ar -dimensional space of states of 
the model. Apparently, the construction of [2] could be modified to resolve this difficulty. We 
hope to address this question in the future. 

It is reasonable to expect that mathematical structures, similar to those described above 
(namely, the partial differential equations and Painleve type recurrence relations), should mani- 
fest itself in other problems, closely related to the 8V- model with rj = n/3. The most immediate 
candidate is the corresponding "eight-vertex solid-on-solid" (8VSOS) model which belongs to a 
rich variety of algebraic constructions associated with the 8V- model [8] . Recently, Rosengren [9] , 
motivated by considerations of the 3-coloring problem [10], studied precisely this 8VSOS-model 
with r\ = 7r/3 in the case of the domain wall boundary conditions. In Eq.(8.11) of his paper [9] 
he introduced a set of two- variable polynomials (defined recursively), related to the partition 
functions of the 8VSOS model on finite lattices. A detailed inspection of these polynomials 
suggests that they satisfy a partial differential equation, which is completely analogous (though 
not identical) to our equation (|44jl ! This new differential equation is presented in Appendix B. 

Furthermore, we found that some one-variable specialization of Rosengren's polynomials, 
also defined in [9], satisfy a recurrence relation, which is extremely similar to the Painleve IV 
type recurrence relation (|48p of this paper. This new relation is also presented in Appendix B. 
As yet it is not written in a canonical form for Painleve IV type relations, however, we expect 
that it could be brought to such form by a suitable change of variables. We also expect that 
this new recurrence relation for the 8VSOS-model can be connected with the Picard elliptic 
solutions of the Painleve VI following the method of our previous paper [1] . It seems it would 
be extremely interesting to further compare our results with those of ref. [9]. 
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Appendix A. Polynomials P n (x,z), s n (z), s n (z), p n (z) and q n (z). 

In this Appendix we present explicit expressions for the polynomials V n (x,z), s n (z), s n (z), 
Pn(z) and q n (z) for small values of the their index n. 

The two- variable polynomials V n (x,z), represented by Eq. ([34"|) . are defined as solutions of 
the TQ-equation (I38p normalized by (|4ip . These polynomials can be efficiently calculated from 
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the differential equation ([44|) . 

V (x,z) = l, Vi(x,z)=x + 3, V 2 (x,z) = x 2 (l + z) + 5x(l + 3z) + 10, 

V 3 (x, z) = x 3 (l + 3z + 4z 2 ) + 7x 2 (l + 5z + I8z 2 ) + 7x(3 + 19z + 18z 2 ) + 35 + 21*, 

7> 4 (a?, 2) = x 4 (l + 6z + 18z 2 + 30z 3 + 9z 4 ) + 9x 3 (l + 8z + 38z 2 + 152z 3 + 57z 4 ) (A.l) 

+ 18x 2 (2 + 192 + lllz 2 + 217z 3 + 99z 4 ) + 12x(7 + 72z + 171z 2 + 198z 3 ) 

+ 18(7 + 14z + llz 2 ). 

The polynomials s n (z) and s n (z), where n E Z, are defined by Eq. (|4"7|) and the recurrence 
relation (|48p . (|49p . For non- negative n they coincide with the coefficients of the highest and 
lowest powers of x in V n (x, z), corresponding to k = n and k = in the expansion (|34p . 

s_ 5 = sfe (81 + 1215 z + 10206 z 2 +64638 z 3 + 353565 z 4 + 544563 z 5 +352836 z 6 ), 

s_ 4 = ^ (27 + 270 z + 1620 z 2 + 7938 z 3 + 3969 z 4 ), 

s_ 3 (z) = ^(9 + 54z + 225z 2 ), 

s- 2 (z) = i(3 + 9z), 

8_i(^)=a («)=aiW = l, (A.2) 

s 2 (z) = 1 + z, 

s 3 (z) = l + 3z + 4z 2 , 

s 4 ( z ) = l + 6 z + 18 z 2 + 30 z 3 + 9 z 4 , 

s 5 (z) = l + 10z + 51z 2 + 168z 3 + 355z 4 + 318z 5 + 121z 6 . 

and 

s (z) = 1, si(z) = 3, s 2 (z) = 10, 
s 3 (z) =35 + 21z, 

a 4 (z) = 126 + 252 z + 198 z 2 , ( A 3) 

s 5 (z) = 462 + 1980 z + 3960 z 2 + 4004 z 3 + 858 z 4 , 

s 6 (z) = 1716 + 12870 z + 47190 z 2 + 105820 z 3 + 143520 z 4 + 90558 z 5 + 24310 z 6 . 
The polynomials p n (y) and q n (y) are defined by the factorization relations (I56p and (I58p in 
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Conjecture E, 

p_ 3 (y) = 1 - 3y + 12y 2 - 30y 3 + 81y 4 - 63y 5 + 66y 6 , 
p- 2 (y) = l-2y + 5y 2 , 

p-i(y) =po(y) = l, 

pi(y) = l + y + 2y 2 , (A.4) 
P2(y) = l + 2y + 7y 2 + 10y 3 + 21y 4 + 12y 5 + lly 6 , 
pa(y) = l + 3y + 15y 2 + 35y 3 + 105y 4 + 195y 5 + 435?/ 6 

+ 555 y 7 + 840 y 8 + 710 y 9 + 738 y 10 + 294 y 11 + 170 y 12 . 

and 

g _ 3 (y) = l + 3y 2 + 39y 4 + 21y 6 , 
g_ 2 (y) = l + 3y 2 , 

<?-i(y) = go(y) = l, 

( A -5) 

gi(y) = 1 + 32/ 2 , 

ga(y) = l + 8y 2 + 29y 4 + 26y 6 , 

g3 (y) = 1 + 15y 2 + 112y 4 + 518y 6 + 1257 y 8 + 1547y 10 + 646y 12 . 

Finally, we list polynomials A n (Q from the expressions (|66j) for the alternative spin compo- 
nents ([65]) . 

^(0 = 1, A 2 (C) = 2, A 3 (C) = 7 + C 2 , ^ 4 (C) = 2(3 + C 2 )(7 + C 2 ), 

A(C) = (3 + C 2 )(143 + 99C 2 + 13C 4 + C 6 ), 

MO = 2(26 + 29C 2 + 8C 4 + C 6 )(143 + 99( 2 + 13( 4 + C 6 ), 

A 7 (Q = (26 + 29C 2 + 8C 4 + C 6 ) 

x (8398 + 14433C 2 + 7665C 4 + 2010C 6 + 240C 8 + 21C 10 + C 12 ), 

(A.6) 

Ag(C) = 2(646 + 1547C 2 + 1257C 4 + 518C 6 + 112C 8 + 15C 10 + C 12 ) 

x (8398 + 14433C 2 + 7665C 4 + 2010C 6 + 240C 8 + 21C 10 + C 12 ) 

A 9 (C) = 2(646 + 1547 C 2 + 1257 C 4 + 518 C 6 + 112 C 8 + 15 C 10 + C 12 ) 

x (1411510 + 4598551 C 2 + 5518417 C 4 + 3530124 C 6 + 1331064 C 8 + 
+ 327810 C 10 + 53382 C 12 + 5820 C 14 + 506 C 16 + 31 C 18 + C 20 )- 

Appendix B. Comments on the 8VSOS-model 

Recently, Rosengren [9], motivated by considerations of the 3-coloring problem [10], studied 
the 8VSOS-model with r) = vr/3 in the case of the domain wall boundary conditions. In Eq.(8.11) 
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of his paper [9] he introduced a set of two- variable polynomials, related to the partition functions 
of the 8VSOS model on finite lattices. Here denote these polynomials as P„ (t,s), adding 
the superscript "SOS" to indicate their relevance to the 8VSOS-modeQ. A detailed inspection 
of these polynomials suggests that 

Conjecture 5. The polynomials Pn SOS \t,s), for even values of n = 0,2,4..., are uniquely 
determined (up to a numerical normalization) by the the following partial differential equation 
in the variables x and s, 

{A^ OS \t,s)d! + Bi sos )(t,s)d t + Ci SOS \t,s) +T( sos \t,s)d s } Pi sos \t,s) = 0, (B.l) 

where 

A (sos) ^ s ) = 2t (1 - t) (1 + 2 s - t) (s + 2 s 2 - 2t - st) (2t + st - s) , 
B ( n SOS \t,s) =-4 (2 + s) 2 t 4 -4 (2 + s)(-2s 2 + s 2 n + 3 sn - 5 s - 3 + n) t 3 

+ (40 sn + 8 s 4 n + 40 s 3 n + 60 s 2 n + 8 n - 8 s 4 - 60 s 2 - 8 - 44 s 3 - 36 s) t 2 

-4 (1 + 2 s) s (-s 2 + 3 s 2 n + 5 sn - 3 s + 3 n - 2) t + 4 ns 2 (1 + 2 s) 2 , 
C n SOS) (t,s) = 2 n (2 + s) 2 (1 + n) t 3 - 4 n (2 + s) (1 + sf (1 + n) t 2 

+ns (8 s 3 + 4 s 3 n + 26 s 2 + 15 s 2 n + 18 s + 18 sn + 2 + 5 n) t 

-4n(l + 2s)s(s + sn + l), 
T( 505 )(t,s) =4 (l-s 2 )s(2 + s)(l + 2s)t . 

(B.2) 

A similar differential equation exists for odd values of n. Obviously, the above property is a 
counterpart of the partial differential equation (|44p in the main text of this paper. 
Next, define one- variable polynomial^! 



A SOS \s) = (l + 2sy ii L2J (l + ^) 4 Pi S ° S) (l + 2s,s) (B.3) 
where [x] denotes the integer part of x. We suggest that 



4 Here we use the variables t and s instead of x and £ used in [9]. These variables are similar, but not identical, 
to our variables x and z in (|23|) . In particular, the variable t (which corresponds to x in [9]) is also connected 
with the spectral parameter, while the variable s (denoted as C, in [9]) is related to the elliptic nome q. 

5 These polynomials are simply related to those introduced in [9], 



^( S ) = (l + 2 S ) l ' i] (l + 



4 \i R) {s) 



where pi^\s) are defined by the first unnumbered equation after the Proposition 3.12 in [9]. We thanks 
H.Rosengren for sending us the modified definition (|B.3|I . 
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Conjecture 6. The polynomials pi^ OS \s) satisfy the following recurrence relations 
s(s - l) 2 (s + 2) (2a + l)d 2 s log p { n SOS \s) + 2{s - l)(s 3 - 3s 2 - 6s - l)d s logp { n SOS) (s) 



4(2n + l)(2n + 3) 



pi+ Q i S) ( s )pi-°i S) ( s ) (B.4) 



2 



+ (22n 2 + 35n + 18) s 2 + (46ra 2 + 98n + 42) s + 13n 2 + 29n + 12 = 
with the initial condition po(s) = 1, pi(s) = 1 + 3s. 

The structure of the relation (|B.4h is very similar to that of the recurrence relation (I48p for 
the tau- functions of Painleve VI equation. It should be noted that Eq. (|B.4p is not written in a 
canonical form for such recurrence relations. Nevertheless, we expect that it could be brought 
to such form by a suitable change of variables. We also expect that (|B.4|) can be connected with 
the Picard elliptic solutions of the Painleve VI following the method of our previous paper [1]. 
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